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ABSTRACT 

Compact group galaxies often appear unaffected by their unusually dense environment. Closer 
examination can, however, reveal the subtle, cumulative effects of multiple galaxy interactions. 
Hickson Compact Group (HCG) 59 is an excellent example of this situation. We present a 
photometric study of this group in the optical (H5T), infrared (Spitzer) and X-ray {Chandra) 
regimes aimed at characterizing the star formation and nuclear activity in its constituent galaxies 
and intra-group medium. We associate Hve dwarf galaxies with the group and update the velocity 
dispersion, leading to an increase in the dynamical mass of the group of up to a factor of 10 (to 
2.8 X 10^^ M 0 ), and a subsequent revision of its evolutionary stage. Star formation is proceeding 
at a level consistent with the morphological types of the four main galaxies, of which two are 
star-forming and the other two quiescent. Unlike in some other compact groups, star-forming 
complexes across HCG 59 closely follow mass-radius scaling relations typical of nearby galaxies. 
In contrast, the ancient globular cluster populations in galaxies HCG 59 A and B show intriguing 
irregularities, and two extragalactic H ii regions are found just west of B. We age-date a faint 
stellar stream in the intra-group medium at ~ 1 Gyr to examine recent interactions. We detect a 
likely low-luminosity AGN in HGG 59A by its ~ 10^° erg s“^ X-ray emission; the active nucleus 
rather than star formation can account for the UV-I-IR SED. We discuss the implications of our 
findings in the context of galaxy evolution in dense environments. 

Subject headings: galaxies: clusters: individual (HCG 59) — galaxies: star clusters — galaxies: evolution 
— galaxies: interactions — galaxies: active — galaxies: dwarf — galaxies: fundamental parameters 
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1. Introduction 

Compact galaxy groups populate the high den- 
sity tail of the galaxy number density distribution. 
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The systems catalogued by Hickson (1982, Hick- 
son Compact Groups, or HCGs) exhibit some fea- 
tures, such as dynamical and evolutionary states, 
elliptical fractions and X-ray properties of the 
intra-group medium (IGM) similar to galaxy clus- 
ters, their massive, more populous counterparts. 
In contrast to the well-studied cluster galaxies, 
however, the specific effects of the compact group 
environment on the evolution of its galaxies are 
not yet clear. 

HGGs are defined through criteria of isola- 
tion and surface brightness^ that give rise to self- 
gravitating, dense groupings of a few (typically 
four) main members. Because of their masses, 
these galaxies orbit around the group barycen- 
ter rather sluggishly, with velocity dispersions 
on the order of ucG 250 kms“^ (Tago et al. 
2008; Gox 2000), cf. galaxy cluster dispersions 
of cTciuster 750 kms“^ (Binggeli et al. 1987; 
The & White 1986; Cox 2000). This trait makes 
HCGs valuable laboratories for galaxy evolution: 
the low velocity dispersions force some galax- 
ies into strong, prolonged interactions while oth- 
ers appear undisturbed but are apparently un- 
dergoing enhanced secular evolution. That is to 
say, this latter population is affected by gravita- 
tional interplay with their neighbors, but evolve 
more subtly, without obvious, strong interactions 
(Konstantopoulos et al. 2010). 

Relating the various observational characteris- 
tics of compact groups to those of clusters is im- 
portant for understanding whether they constitute 
their own class, or if they are simply mini-clusters. 
Perhaps more appropriately, structures like com- 
pact groups may be considered plausible building 
blocks of clusters at higher z (e. g. Fujita & Goto 
2004; Rudick et al. 2006). Revealing past inves- 
tigations of HGGs as a class have focussed on 
gas content. Their members are typically defi- 
cient in H I gas when compared to galaxies of 
similar morphological types and masses (e. g. the 
sample of isolated galaxies in Haynes & Giovanelli 
1984). Verdes-Montenegro et al. (2001) proposed 
an evolutionary sequence based on mapping the 
spatial distribution of H i across a large sample of 
HGGs. Johnson et al. (2007, hereafter J07) added 


^ Spf > 3 9 q, i. e. a circular area defined by three galaxy- 
mean-radii about the group is devoid of galaxies of com- 
parable brightness. A group surface brightness of /r < 
26.0 mag defines galaxy density. 


to this investigation by quantifying the gas rich- 
ness of twelve groups with the relation of H i- 
to-dynamical mass, log(MHi)/ log(Mdyn). This 
gave rise to the hypothesis of an alternate, two- 
pronged evolutionary diagram for HCGs, which we 
explored in Konstantopoulos et al. (2010). In one 
path, the galaxies have strong interactions before 
exhausting their cold gas reservoirs for star forma- 
tion, in the other, gas is processed by star forma- 
tion within individual galaxies prior to late-stage 
dry mergers. 

Furthermore, the mid-IR colors of HCG galax- 
ies show an interesting bimodal distribution that 
distinguishes star-forming from quiescent systems. 
Walker et al. (2010) interpret this statistically sig- 
nificant gap as evidence for accelerated galaxy evo- 
lution in the compact group environment. Their 
similar mid-IR color distributions relate HCGs to 
the infall regions of clusters and set them apart 
from any other galaxy sample compared, interact- 
ing or quiescent. This theme was expanded by 
Tzanavaris et al. (2010) who found this gap appar- 
ent also in the distribution of specific star forma- 
tion rates for HCG galaxies. These observations 
together point to compact groups as local exam- 
ples of the plausible building blocks of clusters in 
the early universe. 

In addition, HCGs, which are isolated by selec- 
tion, could potentially help explain the evolution- 
ary history of some field ellipticals. For example, 
Rubin et al. (1990) originally proposed (see also 
Gallagher et al. 2008) that HCG 31 will evolve 
into a single, field elliptical through a wet merger 
(one where gas is still available during the interac- 
tion). ‘Fossil groups’, the probable ultimate fate of 
isolated groupings, were examined by Jones et al. 
(2003), who defined a criterion of diffuse X-ray 
emission in excess of 0.5 x IO^^/i^q^ erg s“^ for 
such a classification. This arises from the process- 
ing of a group’s IGM during a merger (or series 
of mergers), but the low total mass of most local 
compact groups suggests their potential well lacks 
the depth required to heat the IGM to X-ray de- 
tectable levels (Mulchaey & Zabludoff 1998). Us- 
ing multiple mergers as a vehicle toward a fossil 
group end-state maps one path of galaxy evolution 
from the ‘blue cloud’ of star-forming disk galaxies 
to the ‘red sequence’ of quiescent bulge-dominated 
galaxies (Bell et al. 2004). 

Fossil group formation may provide an analogy 
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to cluster centers or sub-clumps where the buildup 
of cD galaxies occurs. If this turns out to be valid, 
the study of compact groups could also help illumi- 
nate morphological transformations in the inner- 
most cores of clusters. Exploring these different 
scenarios may prove fruitful for our understanding 
of galaxy evolution and the buildup of stellar mass 
in the universe. Making meaningful progress in 
this area requires detailed multi- wavelength stud- 
ies in order to map the range of physical processes 
affecting galaxies that are found in these environ- 
ments, determine their histories, and project their 
evolution. 

A consistent treatment of a large sample of 
HCGs is therefore in order. In this work we 
continue the series of Gallagher et al. (2010) and 
Konstantopoulos et al. (2010) and provide a com- 
prehensive, multi-wavelength study of HCG 59. 

We will look at the current state of the group 
through its star formation and nuclear activity; 
investigate its past through the star cluster pop- 
ulations; try to unravel the history of mergers in 
the group; examine its dwarf galaxy system; and 
place it in the context of HGGs in general. 

The core of HCG 59 consists of four giant 
galaxies, a typical number for HGGs in general. 

The group lies at a distance of 60 Mpc, based 
on a recession velocity of vr = 4047 kms“^ 
(Hickson et al. 1992, corrected to the reference 
frame defined by the 3K Microwave Background) 
and Hq = 73 kms“^Mpc“^. Three of the galax- 
ies, A (type Sa), B (EO), and G (Sc), have 
seemingly undisturbed morphologies, and the 
fourth (D, Im) is an unusually large irregular with 
a normal, peaked light profile. The total stel- 
lar mass of the group is = 3.14 x 10^° Mq 

(from the 2MASS A's-band luminosities; Tzanavaris et al. 
2010), while the Hi mass of Mhi = 3.09 x 10® Mq 
is comparable to the value expected for the mor- 
phological types and stellar masses of the member 
galaxies, according to Verdes-Montenegro et al. 

(2001). This is therefore a somewhat gas-rich 
compact group, given that HGGs typically con- 
tain only about a third of the Hi expected. On 
the other hand, the J07 scheme classifies the Hi 
content of the galaxy group as a Type II, i. e. in- 
termediate in gas content, according to its ratio of 
gas-to-dynamical mass of log(MHi)/ log(Mdyn) = 

0.81 ± 0.05. These classihcations are based on dif- 
ferent criteria and the disparity can thus be rec- 


onciled. Table 1 summarizes some of the general 
characteristics of the four galaxies, while Table ?? 
presents some derived and literature values of the 
mass content and nuclear identifications in the 
four galaxies. 

This paper is organized in the following way: 
Section 2 presents the optical, IR, and X-ray 
datasets used throughout this work. Section 3 pro- 
vides a full account of the young and old star clus- 
ter populations, which we use as our prime diag- 
nostics of current star formation and ancient inter- 
actions. In Section 4 we discuss the main findings 
of this work. Einally, in Section 5 we summarize 
the work presented and offer ties to previous and 
future work in this series. 

2. Observations 

2.1. HST optical imaging 

The analysis presented in this paper is based 
largely on H5T-AGS/WFG multi-band data. Im- 
ages were taken in the F435W, F606W and 
F8I4 W bands in two pointings to cover all known 
giant group members. We will refer to these filters 
as H 435 , Veo6, hi4 (and the set as BVI) to denote 
the closest matches in the Johnson photometric 
system. The notation does not, however, imply 
a conversion between the two systems. The ob- 
servations were executed on 2007 February 24, as 
part of GO program 10787 (PI: Jane Charlton). 
The exposure times were 1710, 1230 and 1065 sec- 
onds in the BVI bands respectively. Three equal 
sub-exposures were taken with each filter with 
a three-point dither pattern (sub-pixel dither- 
ing). Images were reduced ‘on the fly’ to pro- 
duce combined, geometrically corrected, cosmic- 
ray cleaned images. For the analysis of point 
sources, we used the standard HST pipeline prod- 
ucts with a nominal pixel scale of 0'.'05 per pixel. 
For analysis of the extended sources, we ran 
MultiDrizzle (Fruchter & Sosey 2009) with the 
pixel scale set to 0'.'03 per pixel to improve the 
spatial resolution. The absolute image astrometry 
was checked with the world coordinate system of 
the Two Micron All-Sky Survey catalog (2MASS; 
Skrutskie et al. 2006) by identifying four unsatu- 
rated point sources in common; the average offset 
was ^ 0.01" in RA and Dec. The four main galaxy 
/ 814 -band light profiles were fit with Sersic prohles 
using GALFIT (Sersic 1968; Peng et al. 2010a); 
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Table 1: Basic information on HCG 59 main members 


Identifier 

Coordinates®' 

(J2000) 

Type 

H89'= 

m 

(mag) 

VR 

(kms“^) 

References^ 

A: IC 0737 

11:48:27.55 +12:43:38.7 

Sa 

14.82 (B) 

4109 

[1], [2], [3] 

B: IC 0736 

11:48:20.08 +12:42:59.5 

EO'^ 

15.60 (B) 

4004 

[4], [2], [5] 

C: KUG 1145+129 

11:48:32.44 +12:42:19.5 

Sc 

15.90 (g) 

4394 

[4], [2] 

D: KUG 1145+130 

11:48:30.64 +12:43:47.8 

Inr 

16.00 (5) 

3635 

[4], [2] 


‘^Coordinates are the centroids from fitting the HST / 8 i 4 -band images of each galaxy with Sersic profiles. See § 2.1 for more 
details. 

^[1]: Evans et al. (2010); [2]: de Vaucouleurs et al. (1991); [3]: Hickson et al. (1992); [4]: York et al. (2000); [5]: Falco et al. (1999). 
^Hickson et al. (1989) 

‘^The RC3 designation for 59B is‘ SO?’, and so there is some uncertainty as to its classification. 



Fig. 1.— HST BVI color-composite imaging of HCG 59. The red colors of galaxies A and B imply little 
star formation, while C and D show some signatures of nebular emission in green. Also visible is the newly 
catalogued dwarf galaxy I (Section 4.1). 


4 



Table 2: Masses and star formation rates of HCG 59 galaxies 


ID 

(xl09 Mq) 

SFR" 

(Mq yr-1) 

sSFR** 

(xl0-“ yr-i) 

Nucleus ® 

A 

17.40 

10.2 

4.99 ±0.67* 

28.66 

G/AGN 

B 

8.29 

< 8.8 

0.02 ±0.01 

0.19 

G 

G 

3.03 

< 7.2 

0.16 ±0.03 

5.15 

Hii 

D 

2.67 

< 6.6 

0.48 ±0.04 

18.15 

Hii 


masses, star formation rates (SFR) and specific SFRs (sSFR) are drawn from Tzanavaris et al. (2010). The published 
stellar masses were off by a factor of 7.4; these values have been corrected for this error. 

*'From the CO observations of Verdes-Montenegro et al. (1998) 

®Prom Martinez et al. (2010); ‘C’ stands for ‘composite’, i. e. one that falls at the Hll/AGN overlap region, as defined in 
Kewley et al. (2006); see Sections 4.4, 4.5 for discussion on these designations. 

^This value is heavily affected by the AGN in galaxy A, as will be elaborated in Section 2.6. 


the best-fitting centroid positions are given in Ta- 
ble 1. 

We used the images, presented in Figure 1, to 
characterize the optical morphology of the galax- 
ies, and to detect and photometer star clusters 
and cluster complexes. All reported magnitudes 
are in the Vega magnitude system. In Section 3, 
we present the analysis of these two scales (clus- 
ters and complexes) of the star formation hierar- 
chy and also distinguish between young massive 
clusters (YMCs) and globular clusters (GCs). 

2.2. Optical point source photometry 

We follow the same rationale applied in our pre- 
vious work (Gallagher et al. 2010; Konstantopoulos et al. 
2010) and use star clusters to infer the star forma- 
tion activity and history in each of the HGG 59 
galaxies. At the adopted distance to HGG 59 of 
60 Mpc, we expect some contamination by super- 
giant stars, which can have absolute V-band mag- 
nitudes as bright as —8.5 (Efremov et al. 1986). 

At this distance, one AGS pixel measures ~ 13 pc, 

(c/. the average star cluster radius of ^ 4 pc; e. g. 
Scheepmaker et al. 2007), meaning that clusters 
are at most marginally resolved and can be consid- 
ered point sources for the purposes of selection and 
photometry. We select clusters using the method 
described in Gallagher et al. (2010); in brief, we 
perform the initial selection on median-divided im- 
ages, require selection in all three bands, and filter 
the resulting catalog using point spread function 
(PSF) photometry. Our PSF filtering applied the 
following criteria from the output of the ALLSTAR 


routine in IRAF^: y values below 3.0; a sharp- 
ness in the range [—2.0, 2.0]; and a photometric 
error less than 0.3 mag. Aperture corrections are 
first measured between 3 and 10 pixels and then 
added to the Sirianni et al. (2005) corrections to 
infinity. Finally, foreground (Galactic) extinction 
with E{B — V) = 0.037 is accounted for using the 
standard Galactic extinction law (a correction of 
Ay 0.12 mag; Schlegel et al. 1998). 

In order to fortify the selection against stars, 
we apply a conservative absolute magnitude cut 
at My < —9 mag, which produces the high- 
confidence sample. We do, however, define a larger 
sample by relaxing the magnitude cut and apply- 
ing stricter PSF-fitting criteria to detect globular- 
clusters, which are expected to be fainter and 
point-like. In order to minimize contamination 
from marginally resolved sources such as com- 
pact background galaxies, we follow Rejkuba et al. 
(2005) and apply hyperbolic filters (starting nar- 
row for bright sources and widening for fainter 
sources) with a maximum cut-off at the above 
mentioned criteria of magnitude error, y and 
sharpness. We will refer to this as the extended 
sample. 

The application of these criteria assigns 240 
bright star cluster candidates (SGGs) to the high- 
confidence sample and 948 to the extended sample. 
Specifically, the numbers of detected SSGs (ex- 


^ IRAF is distributed by the National Optical Astronomy 
Observatories, which are operated by the Association of 
Universities for Research in Astronomy, Inc., under coop- 
erative agreement with the National Science Foundation. 
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tended sample numbers in parentheses) in galax- 
ies A through D are 7 (29), 77 (213), 13 (63) and 
65 (217), with a further 78 (426) objects coinci- 
dent with what would be the intra-group medium. 
We will provide a full analysis of these cluster pop- 
ulations in Section 3.1. 

In order to test the completeness of the final list 
of SCCs, we used ADDSTAR to add 3000 artihcial 
stars to the image (over the entire field, includ- 
ing the galaxies) in the apparent magnitude range 
24-28, i. e. absolute magnitudes of (—9.89, —5.89). 
Because the final catalogue only contains sources 
detected in all three filters, we include this effect 
by calculating completeness fractions based only 
on artihcial stars detected in all three bands (e. g. 
Da Rocha et al. 2002). The limiting magnitudes 
for the 90% and 50% recovery rates are (26.56, 
27.25), (26.51, 27.19) and (26.47, 27.15) in the 
B 435 , Veoe and Jgi 4 bands respectively (after pho- 
tometric corrections are applied) . For the distance 
modulus used of 33.89 mag, m = 26.5 mag corre- 
sponds to M ~ —7.4 mag. 

Our assessment of the state of star formation 
in HCG 59 is not limited to star clusters. Star 
cluster complexes represent a larger scale of star 
formation, as the optically blended concentrations 
of gas, stars, and dust that make up small star- 
forming regions, and likely include groups of clus- 
ters. In contrast to star clusters, these can be re- 
solved to even greater distances than studied here, 
as the fractal distribution of gas about a galaxy 
gives rise to such structures at all scales (e. g. as 
demonstrated for M33 by Bastian et al. 2007). 

2.3. Globular Cluster Candidate Selection 

Globular clusters are also selected from the 
HST images. Since the process is tuned to the 
color distributions found in HGG 59, we provide a 
full account below. As contamination from super- 
giants is not a problem for objects with GG-like 
colors, we adopt a fainter magnitude limit to se- 
lect old GG candidates than we used for SGGs. 
We have chosen a cutoff at Veoe = 26, which cor- 
responds to My —7.7 at our adopted distance 
modulus for HGG 59, or slightly more luminous 
than the expected peak in the globular cluster lu- 
minosity function (GGLF) at My ~ —7.4 (e.g. 
Ashman & Zepf 1998; Harris 2001). The major- 
ity of GG candidates brighter than this limit lie 
above the 90% photometric completeness level in 


all 3 filters. Assuming a Gaussian GG luminosity 
function with a peak at My = —7.4 ± 0.2 and a 
dispersion a = 1.2 ± 0.2, our faint-end cutoff then 
samples 39 ± 8% of the entire GGLF. 

We have selected GG candidates (GGGs) ac- 
cording to the color-space distribution of Milky 
Way GGs. We de-reddened the colors of globular- 
clusters from the Harris (1996) catalog by their 
listed E{B — V) values, and then defined a paral- 
lelograiu based on the intrinsic {B — V) — {V — I) 
color distribution of the MW GGs. This parallelo- 
gram was then converted to the AGS filters using 
the ‘synthetic’ transformations in Sirianni et al. 
(2005). The color selection region adopted here 
is 0.10 mag wider in {V — I) than that used in the 
analysis of HGG 7 (Konstantopoulos et al. 2010), 
but still does not exceed the boundaries of the 
Harris (1996) MW GGs. To quantify, 95 of 97 
MW GGs from the Harris (1996) catalog (those 
with BVI information) lie within this box. All 
point sources in the HGG 59 fields with 1 a er- 
ror bars that overlap our color selection region are 
considered GG candidates, and are plotted in Fig- 
ure 2. 

Due to the close (projected) proximity of the 
galaxies in the group, it is likely that the halo 
GGs in each system will appear superposed. In an 
attempt to quantify the GGGs in each galaxy, we 
use the relationship between the galactic mass and 
the radial extent of the GG systems in galaxies of 
Rhode et al. (2007). To compute the expected size 
of each halo, we have adopted the mass-to-light 
conversions in that work, although we stress the 
general conclusions we reach are not dependent 
on the detailed size of any given halo. As the pre- 
dicted masses of all of the group galaxies are just 
below the lowest mass galaxies in the Rhode et al. 
(2007) sample, we adopted a radial extent of 15 
kpc (or 56" at the assumed distance to HGG 59) 
for each of galaxies A, B and G. Galaxy D is a 
lower luminosity system, and as seen in Figure 2, 
GGGs in this object already lie within the pro- 
jected halo of GGGs in galaxy A. Discussion of 
the individual GG systems in each group galaxy 
follows in Section 3.3. 

2.4. Background and Foreground Contam- 
ination 

Gontamination in our color-selected sample 
of GG candidates is expected from a variety of 
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Fig. 2. — Positions of detected globular cluster candidates, marked on HST I^u imaging. We have also 
marked the expected GC haloes with blue boundaries. ‘Outer B’ refers to the location of a possible excess of 
GCs, discussed in the text. This is symmetric to the ‘A-B bridge’ (with respect to B; see Fig. 4) and could 
be a tidally induced redistribution of the GG population (see Sections 3.3 and 4.4). 
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sources, including foreground Milky Way halo 
stars, reddened young clusters, and unresolved 
background galaxies. With the small number of 
GC candidates present in some of the HCG 59 
galaxies (discussed below), contamination can be 
significant. 

Predictions from Milky Way star count mod- 
els (the Besangon model of Robin et al. 2003) sug- 
gest that only 3-4 foreground Milky Way stars will 
appear in the magnitude and color ranges for ex- 
pected GCs in each of our ACS fields. Determining 
the contamination from younger, reddened clus- 
ters is more difficult, particularly in the central re- 
gions of the late- type galaxies C and D, where such 
objects could be present. Unresolved background 
galaxies are not likely contributing in any signifi- 
cant way to the numbers of objects in our helds; 
analyses of the background objects (with GC-like 
colors) in HCG 7 (Konstantopoulos et al. 2010) 
showed that the predicted foreground Milky Way 
star counts were similar to the observed putative 
background contamination, leaving little room for 
background galaxies to contribute significantly. 

We also consider the Pirzkal et al. (2005) anal- 
ysis of stars in the Hubble Ultra-Deep Field 
(HUDF). Within the range of colors shown in 
Figure 8, they found the main contaminant of 
‘void sky’ to be M-stars, however, with a Veoe- 
/si 4 of ~ 2.0, they are too red to be considered 
in our analysis. All Main Sequence stars detected 
in the HUDF are too bright to be mistaken for 
star clusters by our detection algorithm. In fact, 
the only class of stellar object that can be found 
in the color-space occupied by our cluster can- 
didates is white dwarfs, which Pirzkal et al. find 
to have a density of 1.1 ± 0.3 x 10“^ pc“^. The 
maximum Galactic volume covered by our two 
pointings is a cube of ^ 8000 px on a side, or 
^ 0.24 pc“^, assuming a scale height of 400 pc (the 
maximum height quoted by Pirzkal et al.). Such 
a volume might be expected to host ^ 2.6 x 10“^ 
white dwarfs. Therefore, we consider this poten- 
tial source of contamination to be negligible. 

The dominant source of contamination to our 
GCG sample will actually be stars from the Sagit- 
tarius dwarf galaxy. Our ACS fields (at I 254°, 
b ^ -1-69°) are superposed on a rather dense part 
of the tidally induced (and bifurcated) leading 
arm (stream ‘A’ from Belokurov et al. 2006) of 
Sgr (e.g., Majewski et al. 2003; Belokurov et al. 


2006; Newberg et al. 2007; Yanny et al. 2009); 
the presence of such streams is not accounted 
for in traditional Milky Way star count mod- 
els. To investigate the impact of Sgr Stream 
stellar populations, in Figure 3 we have overlaid 
12 Gyr isochrones of Marigo et al. (2008) with a 
range of metallicities expected for the Sgr lead- 
ing arm, [M/H] ^ — 1 ± 0.5 (e.g. Chou et al. 
2007; Yanny et al. 2009), at distances between 
26 and 36 kpc onto color-magnitude diagrams 
of the point sources in both of our ACS fields 
(assuming d ~ 31 kpc, with a spread of ±5 
kpc; Newberg et al. 2007; Niederste-Ostholt et al. 
2010; Correnti et al. 2010). From this, we see that 
stars just below the Sgr Stream main sequence 
turnoff do have colors and magnitudes similar to 
that of the brighter (Veoe < 25) GC candidates in 
our study, making some contamination likely. 

To estimate the total contamination in our GCG 
sample, we assume those GC candidates that lie 
far outside the GC system halos (as shown in Fig- 
ure 2) are instead contaminating sources. The 
one exception to this is a region (called ‘outer B’ 
in Figure 2) that lies outside the GC system of 
galaxy B, opposite to the direction of galaxy A. 
We will return to this feature below. There are a 
total of 18 objects in 8.4 arcmin^, or a background 
surface density Ehacfe = 2.3 ± 0.5 arcmin”^. This 
is much higher than the predicted surface den- 
sity of MW halo stars from the Besangon model 
(Emw ~ 0-4 arcmin“^), indicating that Sgr lead- 
ing arm stars are the dominant foreground source 
of contamination in our sample. 

Of course, for the above analysis we are mak- 
ing the assumption that these contaminating ob- 
jects are not bona-fide ‘intra-group’ GCs that lie 
far outside the main galaxies of the group. To test 
this, we compare the Vqoq luminosity function for 
the background source sample with the luminos- 
ity function of the large GC candidate population 
surrounding galaxy B. We show this in the right 
panel of Figure 3. The extrahalo luminosity func- 
tion does not show a sharp rise with increasing 
magnitude as expected of a GC luminosity func- 
tion and exhibited by the GCCs in galaxy B. Al- 
though a definitive comparison is not possible with 
so few objects in the IGM area, this is consistent 
with these IGM objects being contaminants and 
not a part of a diffuse population of intragroup 
GCs. Thus we adopt the surface density above 


HCG 59 



(^435 ^603)0 (^606 ^814)0 



Fig. 3. — Diagnostic diagrams related to the GC population. On the left we show color-magnitude diagrams 
of all detected GCCs. The solid and dashed lines show Marigo et al. (2008) evolutionary tracks for 12 Gyr- 
old main sequence stars, with Z = —0.6 and —1.3 respectively. The lines are double to bracket the distances 
appropriate for the Sagittarius dwarf, with its leading spiral arm superposed along this line of sight. Their 
distribution in the GMD dispels our worries of significant contamination. The right panel shows the lumi- 
nosity function of globular clusters candidates (filled circles) in galaxy B and sources considered ‘background’ 
(open circles and dashed lines). The clear discrepancy supports the quality of our GC selection. The LF of 
‘outer B’ sources is designated with a solid line. 
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as that of the ‘background’ in the analyses that 
follow. 

2.5. Las Campanas wide-field imaging: 
low surface brightness light and dwarf 
galaxies 

We extend the coverage of the HST ob- 
servations through wide-held imaging with the 
Las Campanas Observatory (LCO) 2.5-meter tele- 
scope. We took B- and i?-band images of a 25' 
diameter around the group with the Wide Field 
Reimaging CCD Camera (WFCCD). The data 
were obtained on 2007 July 07 as part of an imag- 
ing campaign that covers all 12 HCGs in the J07 
sample. The B and R hlter exposure times were 
300 s and 600 s, respectively. 

These images allow for the detection of low sur- 
face brightness features, such as the signatures of 
past interactions, over a large area. We present 
this analysis in Figure 4, where we stacked the B 
and R images and applied a Gaussian smoothing 
hlter to the result. This image shows only fea- 
tures and i?-band contours that register at least 
3 (T above the background. We hnd two faint fea- 
tures, a ‘bridge’ that appears to connect galaxies 
A and B and an arc extending from B toward a 
compact structure to its north-west (we will later 
refer to this as the ‘B-I arc’). There is another 
compact, extended source in the space between 
galaxies C and D. 

The original purpose of the LGO observing 
program was to prepare a sample of dwarf galaxy 
candidates for spectroscopic follow-up. Though 
our redshift survey has yet to cover HCG 59, it is 
covered by the Sloan Digital Sky Survey (SDSS; 
York et al. 2000): a spectroscopic search sweep- 
ing a radius of 30 arcminutes around the nom- 
inal center of the group (the geometric center 
of the region enveloping the four known mem- 
bers) yields seven spectra with redshifts in the 
range 0.01-0.02: galaxies C and D and five com- 
pact galaxies. We therefore consider the mem- 
bership of SDSS galaxies J114817. 89-1-124333.1 
and J114813. 50-1-123919. 2, which are covered by 
our wide-held imaging and J114930. 72-1-124037. 5, 
J114940.11-kl22338.6 and J114912. 21-^123753.8 
which lie at projected distances greater than 13 ar- 
cminutes from the group center. The hrst of these 
galaxies is also present in the HST imaging, but 
lies partly in the ACS chip gap. In Section 4.1 we 


attempt to determine whether these are HCG 59 
members through a phase-space analysis. 

The Hickson (1982) naming convention assigns 
letters in order of brightness. Since our imaging 
does not cover all hve dwarf candidates, we used 
the SDSS r-band photometry to consistently clas- 
sify the galaxies as HCG 59 F through J. We have 
omitted the letter E, as it was assigned in the 
original catalog to a background galaxy. We at- 
tempted to measure stellar masses for these galax- 
ies using 2MASS AT-band images (Skrutskie et al. 
2006), however, they are below the detection limit 
of that survey. Table 3 summarizes all of the in- 
formation presented in this section: measured and 
SDSS photometry, radial velocities, galaxy mor- 
phologies and projected distances from the group 
barycenter. The latter two properties will be dis- 
cussed in Section 4.1 

2.6. Spitzer observations: infrared spec- 

tral energy distributions 

The optical imaging was complemented by 
Spitzer imaging in the mid-infrared (IRAC 3.6- 
8 pan and MIPS 24 pan observations) presented 
in JO 7 and shown in Figure 5. In addition to 
the Rayleigh- Jeans tail of stellar photospheric 
emission, the IRAC bands probe the presence 
of hot dust and polycyclic aromatic hydrocar- 
bons (PAHs), while the 24 pan observations trace 
cooler thermal dust emission. The dust and PAH 
emission are both stimulated by star formation 
activity. The harder spectra of active galactic nu- 
clei typically destroy PAH molecules while heat- 
ing dust to hotter temperatures than found in 
galaxies with star formation alone. At low AGN 
luminosities, the IR SEDs are often ambiguous 
(particularly in the presence of star formation; 
e.g., Gallagher et al. 2008). 

The Spitzer images were combined with JHKs 
observations from 2MASS (Skrutskie et al. 2006) 
to plot the IR spectral energy distribution (SED) 
of each galaxy (following J07), presented in the 
frequency-space plot of Figure 6. We have used the 
Silva et al. (1998) templates for galaxies of vari- 
ous morphological types. These map the SED of 
different galaxies as the sum of starlight and gas 
and dust emission from star formation and inter- 
stellar cirrus. We calculate the spectral index of 
the SED within the IRAC bands through a simple 
power-law fit. This was defined by Gallagher et al. 
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Fig. 4. — Co-added, smoothed, background-suppressed B R image of HCG 59 from the LCO duPont 
telescope. The image is aligned to the world coordinate system. The contours mark i?-band fluxes between 
3 and 10 cr above the background (as marked by the circled numbers), while the smoothed B -\- R flux shows 
material brighter than 3 cr above the background level. The asterisk label denotes the foreground star to the 
east of galaxy D, while the double asterisk marks the Hii regions west of B. The horizontal banding north 
of the main group is an artifact in the i?-band image. We detect a streak of material apparently connecting 
galaxies A and B in our field of view. We argue in Section 4.3 that this is a bona-fide tidal bridge, although 
its low surface brightness prohibits an in-depth photometric study. A faint arc of light to the west of galaxy B 
- containing the H ii regions mentioned above (see § 3.1) - appears to connect it with the dwarf galaxy I. 
We discuss the importance of these features in Section 4.3. 
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Table 3: Dwarf galaxies in HCG 59 


ID 

SDSS ID 

(Plate, MJD, Fiber) 

Morphology 


(mag) 

Mr^ 

b 

vr 

(kms“^) 

dec 

(kpc) 

F 

J114930. 72+124037.5 

(1609, 53142, 464) 

Sd 

15.83 

0.09 

-18.15 

3984 

290 

G 

J114940. 11+122338.6 

(1609, 53142, 465) 

dir 

17.34 

0.08 

-16.63 

3195 

471 

H 

J114912.21+123753.8 

(1608, 53138, 625) 

dir 

17.39 

0.08 

-16.58 

4011 

225 

I 

J114817.89+124333.1 

(1608, 53138, 586) 

dim 

17.43 

0.10 

-16.56 

4038 

32 

J 

J114813.50+123919.2 

(1608, 53138, 593) 

dE 

17.69 

0.10 

-16.30 

3942 

83 


“SDSS model r magnitudes. These values have not been corrected for Galactic extinction. 

^’Galactic extinction and recession velocities are from SDSS DR7. 

‘^Absolute r magnitudes assuming a distance modulus of 33.89 and correcting for Galactic extinction. 




Fig. 5. — Spitzer 3.6-8 /tm images of HCG 59. We present these individually to contrast the features evident 
in different bands and in a combined image to view the entire mid-IR picture. Galaxy A is extremely bright in 
the hot dust/PAH-dominated 5.8 and 8 fim bands. Galaxy B is faint in those bands, owing to little emission 
from heated dust. C is consistently bright in the four bands, given its ongoing star formation. D appears 
clumpy at long wavelengths, indicating patchy dust, heated by forming stars. In the 3.5 and 4.6/im images, 
we also see faint, diffuse emission in the region between A and B, which we study in Sections 2.5 and 4.3. 
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(2008) as aiRAC, and it serves as a measure of star 
formation activity. In logarithmic frequency units, 
the flux difference from 8 to 4.5 fim leads to a posi- 
tive gradient in quiescent environments, while star 
formation registers as a negative slope. For HCGs, 
the steepness of the slope is sensitive to the specihc 
star formation rate (Tzanavaris et al. 2010). 

In brief, the SEDs of galaxies B and C follow 
their morphological types of E/SO and Sc, while 
the irregular nature of D does not allow for a 
template to be assigned (we use Sc in Figure 6). 
Galaxy A, nominally an Sa, shows strong excess 
light at wavelengths associated with PAH and/or 
hot dust emission (A > 5.8 ^m). If the flux in 
this region is associated with star formation, an 
Sc template would be more appropriate. This is, 
however, inconsistent with the optical morphol- 
ogy of A. Furthermore, the shape of this emission 
is also consistent with the AGN dust bump (e.g., 
Elvis et al. 1994; Gallagher et al. 2008), so in Sec- 
tion 4.5 we will examine the nuclear activity of 
this galaxy. 

The SED-fitting process could not be car- 
ried out for the new dwarf galaxies catalogued 
in Section 2.5 due to the limited spatial cover- 
age of the Spitzer imaging and the faintness of 
the dwarfs (they were not detected in 2MASS). 

As of July 2011, the Wide- field Infrared Survey 
Explorer (WISE) photometric catalogue did not 
cover HCG 59. 

2.7. The high energy picture: Chandra- 

ACIS observations 

The dataset is completed by Chandra-AClS 
data in the 0. 5-8.0 keV range. Results reported 
here are drawn from the work by Tzanavaris et al. (2011 
in prep.). HGG 59 was observed by Chandra be- 
tween 2008-04-12 and 2008-04-13 at the aim point 
of the back-illuminated S3 GGD of ACIS in very 
faint mode with an exposure time of 39 ks (ob- 
servation ID 9406, sequence number 800743, PI 
S. Gallagher). The data were processed using 
standard Chandra X-ray Center aspect solution 
and grade hltering, from which the level 2 events 
file was produced. Figure 7 shows an adaptively 
smoothed 3-band X-ray image, with optical (/si 4 ) 
contours overplotted for comparison. 

Wavdetect, the CIAO 4.1.2^ wavelet detection 

®http:/ /cxc. harvard.edu/ciao 



13.0 13.5 ' 4.0 ' 4.5 

log(i/) (Hz) 


Fig. 6. — Near-to-mid IR SEDs for the primary 
HGG 59 members. The photometric data, shown 
as diamonds, are drawn from 2MASS {JHKs) and 
Spitzer MIPS/IRAG 3.5-24 ^m (from J07). We 
annotate each panel with the morphological type 
from the literature (top row), as well as the type 
chosen to most accurately represent the data (bot- 
tom row). The bracketed number quotes the age 
in Gyr of the template for the elliptical, or the 
inclination in degrees of the spirals. We also cite 
the value of the aiRAC diagnostic (Gallagher et al. 
2008), which is fit on the wavelength range indi- 
cated by the dashed line. The solid curves repre- 
sent GRASIL templates (Silva et al. 1998) appro- 
priate for the morphological type of each galaxy 
(see text). Galaxies B and G are well represented 
by their nominal templates, while A and D de- 
viate: in the case of galaxy D this is due to its 
irregular nature, while we consider the nuclear ac- 
tivity of A to be the source of the excess mid-IR 
emission (see Section 4.5 for details). 
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Fig. 7. — Adaptively smoothed three-band X-ray image of HCG 59, with overplotted / 814 -band contours. 
Red, green and blue colors represent soft (0.5 — 2 keV), medium (1 — 4 keV) and hard (2-8 keV) X-ray 
emission. Galaxy C (Sc) has no detectable X-ray emission above the background level (the source at the 
limit of its optical isophotes is a likely background galaxy), while D (Im) shows some soft diffuse emission. 
The point source detected in the nuclear region of galaxy A (Sa) is detected in the full band with Ax(o. 5 - 8 . 0 ) = 
1.1 X 10“^° erg s“^. This is surrounded by diffuse emission. We resolve two point sources in the central region 
of Galaxy B (E/SO) with Tx(o. 5 - 8 . 0 ) = (1-4, 1.7) x 10^® erg s“^. Section 4.5 and Figure 17 provide details. 
No diffuse emission is detected from the intra-group medium. 
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tool (Freeman et al. 2002) was used in the soft 
(0.5-2.0 keV), hard (2.0-8.0 keV) and full (0.5- 
8.0 keV) bands to detect candidate point sources. 
The 1024x1024 S3 chip field was searched with 
wavdetect at the 10“® false-probability thresh- 
old. Wavelet scales used were 1, 1.414, 2, 2.828, 
4, 5.657 and 8.0 pixels. Source lists produced by 
wavdetect for each band were matched against 
each other by means of custom-made scripts (K. D. 
Kuntz, priv. comm.) to calculate a unique posi- 
tion for each candidate point source, taking into 
account the varying size of the Chandra PSF 
across the S3 CCD. 

Point source photometry was carried out for 
the objects in the source list using ACIS Extract"^ 
(Broos et al. 2010) which takes into account the 
varying Chandra PSF accross the CCD. Pois- 
son ±1 (T errors on net counts were calculated by 
means of the approximations of Gehrels (1986). 
Sources with measured net counts smaller than 
the 2 cr error were flagged as non-detections. Note 
that this method produces very similar results to 
choosing a binomial probability threshold of 0.004 
in ACIS Extract (Xue et al. 2011). 

As sources have too few counts for reliable spec- 
tral fitting, we apply the method of Gallagher et al. 
(2005) to obtain a rough estimate of the spectral 
shape by using hardness ratios, defined as HR 
= {H — S)/{H + S), where H and S represent 
the counts in the hard and soft bands, respec- 
tively. Briefly, we use the X-ray spectral modeling 
tool XSPEC (Arnaud 1996), version 12.5.0, to sim- 
ulate the instrumental response and transform ob- 
served HR values into an effective power law index 
F (where fx cx E~^ photons cm“^ s“^ keV“^), 
and also obtain associated X-ray fluxes and lu- 
minosities. This modeling includes neutral ab- 
sorption from the Galactic A^h of 2.6 x 10^° cm“^ 
(Kalberla et al. 2005). 

We estimate flux limits of fx ^ 2.7 x 10“^® 
erg cm“^ s“^ (0.5 - 2.0 keV) and fx ^ 1.6 x 
10“^® erg cm“^ s“^ (2.0 - 8.0 keV), corresponding 
to luminosity limits of Lx = 1.3 x 10®® erg s“® 
(0.5 - 2.0 keV) and Lx = 7.7 x 10®® erg s"® 
(2.0 - 8.0 keV). Assuming these limits, we use the 
log N — log S relation of Cappelluti et al. (2007) 
to estimate the number of background sources we 
would expect to detect in the HCG 59 field over 


^http: / /www. astro. psu.edu/xray/docs/TARA/AE.html 


the 1.96 X 10“^ square degree area of the S3 chip. 
This number is 25 and 20 in the soft and hard 
bands, respectively, with a ^ 20% uncertainty. In 
the much smaller area (3.66 x 10“^ square degrees) 
covered by our galaxies, we expect < 1 background 
source in each band. 

We detect a total of 40 sources in the soft band 
and 33 sources in the hard band over the ACIS S3 
field. We thus expect about 15 soft and 13 hard 
sources to be point sources associated with HCG 
59. We find that 9 soft and 3 hard sources are 
located inside the boundaries of the MIR-based 
HCG 59 galaxy regions of J07. We note that 11 
sources that have only hard-band emission are lo- 
cated far from the HCG 59 galaxies and are thus 
likely background AGN. 

In the central regions of two group galaxies, 
there are three notable X-ray point sources de- 
tected with high significance. As these point 
sources are all with 5' of the Chandra optical axis, 
the X-ray positions are the ACIS Extract “mean 
positions of events within the extraction regions” 
® with intrinsic positional uncertainties of a few 
tenths of an arcsecond. There is some additional 
uncertainty from matching the absolute reference 
frames of Chandra and HST, but this is expected 
to be small as both are consistent with 2MASS at 
the 0.1" level. 

The first source, in galaxy A, has a full-band 
luminosity of Tjc(o. 5 - 8 . 0 ) = 1-1 10^® erg s“® and 

an estimated F = 1.3 ± 0.3. The X-ray posi- 
tion of this source is 0.7" from our quoted opti- 
cal position (Table 1). This isolated point source 
in the nuclear region of galaxy A has a luminos- 
ity that is consistent with known low-luminosity 
AGN (e.g.. Ho et al. 2001) and significantly higher 
than individual, luminous X-ray binaries. The 
two point sources found in the central region of 
galaxy B have X-ray positions 0.2" and 1.3" from 
the HST Tband galaxy centroid position (Ta- 
ble 1). Unfortunately, these sources have fewer 
than 10 counts in each band, precluding even a 
rough F estimate. Their full-band luminosities are 
Ax(o. 5 - 8 . 0 ) = (1.7, 1.4) X 10®® ergs”®, respectively. 

The Chandra images are also sensitive to diffuse 
emission from MK degree gas. As can be seen in 
our adaptively smoothed image, some soft, diffuse 


®See Section 5.3 of the ACIS Extract User Manual - 
http: / /www2. astro. psu.edu/xray/docs/TARA/ae_users_guide.html 
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emission is detected in galaxy D (Im), likely asso- 
ciated with star formation, as well as in galaxy A. 
In both cases, the diffuse emission covers an area 
several times the size of the Chandra PSF at that 
location. We obtain an upper limit on the IGM 
surface brightness as follows. We calculate the 
count rate in a source-free region between group 
galaxies, and estimate the corresponding flux for 
kT = 0.5 keV thermal emission using PIMMS.® 
We thus estimate the IGM surface brightness to 
be ^ 7.3 X 10“^^ erg cm“^ s“^ arcsec“^. Finally, 
we note that dwarf galaxy I is not coincident with 
any detected X-ray sources; none of the other new 
dwarf members of HCG 59 are within the Chandra 
field of view. 

We will discuss the implications of these obser- 
vations further in Section 4. 

3. The young and old star cluster popu- 
lations: star formation over a Hubble 
time 

3.1. Young star clusters and the past ^ 
Gyr of star formation 

The population of star clusters is representa- 
tive of star formation as a whole in any system 
(e.g. Bressert et al. 2010). They are formed en 
masse after large events (e.g., Tranche et al. 2007; 
Konstantopoulos et al. 2008, 2009; Konstantopoulos 
2009; Bastian et al. 2009) and at a slower pace 
at all times when a galaxy is forming stars 
(Lada & Lada 2003). The extreme brightness of 
young clusters makes them detectable to large 
distances and therefore a reliable tracer of the 
star formation history of their host galaxy over a 
Gyr or so. Beyond that point in time they are 
referred to as intermediate age clusters and even- 
tually globular clusters. As a whole, the cluster 
population of a galaxy can reveal its star-forming 
history over a Hubble time. 

In this section, we analyze the cluster popu- 
lations of the four main galaxies in HCG 59 us- 
ing the HST PSF photometry described in Sec- 
tion 2.2. We use color-color (CC) and color- 
magnitude (CMD) diagrams to roughly age-date 
the clusters by comparing them to evolutionary 
tracks. The BVI filter combination of our HST 
images lacks coverage below 4000 A, which is 


®http: / /heasarc. nasa.gov/Tools/w3pimms.html 


crucial to breaking the age-reddening degeneracy 
(owing to the inclusion of the Balmer jump and 
near-UV continuum). However, it is still possible 
to infer the passing of intense bursts of star for- 
mation via the clumping of data-points along the 
evolutionary track, and unreddened young clusters 
are clearly evident. 

Figures 8 through 11 show the high-confidence 
sample {i. e. clusters with My < — 9 mag; see Sec- 
tion 2.2) as solid dots, while the extended sam- 
ple is marked with open triangles. The solid and 
dashed red lines (running top to bottom) show 
Marigo et al. (2008) evolutionary tracks for simple 
stellar populations (SSP) of | Zq and Zq. This is 
slightly different from our previous work, where we 
used Bruzual & Gharlot (2003) models. We made 
this change because the Marigo et al. (2008) mod- 
els seem to provide a good fit to both young clus- 
ters and globulars, unlike other model suites which 
focus on one part of the cluster population. The 
green lines that run more or less horizontally show 
Starburst99 (SB99; Leitherer et al. 1999) tracks of 
the same metallicities. These also include nebular 
emission, which we expect to often be present dur- 
ing the first ~ 10 Myr of evolution. At this state 
the cluster is still surrounded by residual gas from 
the time of its formation ionized by UV photons. 
This short-lived phase ends when the first stars 
evolve and explode in supernovae that expel the 
gas. Glusters with colors redder than 0.8 in both 
axes are most likely GGs, although they might be 
highly reddened young clusters. 

When contrasting the extended and high- 
confidence samples, the former appears to spread 
more in color-space. This reflects a mass ef- 
fect intrinsic to star clusters, rather than in- 
dicating contamination. Two recent studies, 
Silva-Villa & Larsen (2011) and Popescu & Hanson 
(2010), independently reached the conclusion that 
lower mass (fainter) clusters, M < 10^ Mq, of- 
ten exhibit deviations from the theoretical model 
tracks. In this mass regime, the underlying stellar 
initial mass functions (IMF) are under-sampled 
and stochastic effects dominate the overall light. 
Given that the IMF is populated randomly, it is 
physically equivalent to creating either one high 
mass star, or ~ 100 low mass stars (given the IMF 
slope). In contrast, a high-mass cluster will popu- 
late the IMF fully. Gonsequently, a population of 
high-mass clusters will have smaller intrinsic pho- 
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tometric dispersion. A cluster with lower mass 
will run out of material before the IMF becomes 
fully sampled, thus leading to the presence of gaps 
and spikes. Because of this effect, the presence of 
a high mass star in a low mass cluster will make 
the cluster appear to have a larger photometri- 
cally derived mass than one with only lower mass 
stars. 

We now treat the population of each galaxy 
individually. In the following paragraphs, all 
star and globular cluster candidates discussed are 
taken from the high-confidence sample. The dis- 
tinctions of young, intermediate, and globular are 
inferred from the locations of candidates along 
evolutionary tracks within the i? 435 -V 606 vs. Vqoq- 
/si 4 color-color space. 

Galaxy A (Sa) hosts a very small detectable 
population of five GCCs and two intermediate- 
age cluster candidates. The CMD shows them 
all to be consistent with having high masses, with 
log(M/M 0 ) > 5. The lack of young clusters, com- 
bined with the high masses of the ones detected, 
indicates that galaxy A has stopped forming stars 
at a high rate. Because the mass-to-light ratio 
of a stellar population increases with age (i e., 
the older a cluster becomes, the higher its mass 
must be for detection), young clusters should dom- 
inate the cluster population when the SFR is high. 
The non-detection of any young clusters is fur- 
ther evidence that the UV-I-IR emission is domi- 
nated by AGN continuum, and the inferred SFR 
of 5 Mq yr“^ is severely overestimated. 

Galaxy B (E/SO) hosts mostly red clusters; 
the GCCs vastly outnumber young cluster candi- 
dates in this system. In fact, the four SCCs in the 
high-confidence sample do not strictly belong to 
galaxy B, but are found in the IGM: two nebular 
sources are located in a clump to the west of the 
galaxy, which we are treating as part of a stream 
that connects this galaxy to the dwarf HCG 591 
in our field of view. We will return to these Hii 
regions in Section 4.3. The two non-nebular SCCs 
are distinct clusters in dwarf galaxy I. This means 
that the presence of patchy dust (as seen in the 
HST images of 59B) by itself is not indicative 
of star formation at a high enough level to pro- 
duce massive clusters, consistent with its SFR of 
0.02 Mq yr“^. 

Galaxy G (Sc) shows a continuous star forma- 
tion history, evident in the smooth distribution of 


datapoints along the evolutionary track. The pres- 
ence of nebular sources indicates some current star 
formation, while the low masses (~ 10^ Mq) de- 
rived throughout the sample imply an overall low 
level of star formation over time. This is in ac- 
cord with the low value of 0.16 Mq yr“^ for the 
SFR of the galaxy. There is no pronounced GC 
population. 

Galaxy D (Im) is an unusually large irreg- 
ular galaxy. It shows a continuous star forma- 
tion history through to the present. The CMD 
shows a handful of SCCs with log(M/MQ) > 5. 
At 0.48 Mq yr“l, this galaxy has the largest SFR 
of those with young clusters. There is no old com- 
ponent in the cluster distribution, no evident halo 
of GCCs, again possibly due to the low mass of the 
system which implies a small globular cluster sys- 
tem. Furthermore, the extended sample does not 
reveal a tight, correlated color distribution charac- 
teristic of GCs. The youngest clusters lie at a typ- 
ical color-space distance of ^ 0.3 mag away from 
the nebular model track along the reddening vec- 
tor, indicating the presence of dusty star-forming 
regions. 

3.2. Star cluster complexes 

One step above star clusters in the star for- 
mation hierarchy is cluster complexes, large ag- 
glomerates of young stars, arranged in a fractal 
distribution that follows the collapse of the pro- 
genitor gas. These structures can be used to un- 
derstand the global star formation activity in a 
galaxy. They are found to be more compact at 
higher redshifts than in the local universe (with 
two to hve times higher mass surface density, as 
found by Elmegreen et al. 2009), although in one 
local interacting compact group, HCG 31, we find 
complexes to be similar to those at intermediate 
redshifts (Gallagher et al. 2010). We argue in this 
series of papers that compact groups might process 
gas more efficiently when interacting than most 
other environments apart from the infall regions 
of galaxy clusters (Walker et al. 2010). HCG 31 
fits well in that context. 

Star-forming complexes are extended amor- 
phous regions, with dust scattered across their 
surfaces. Their boundaries were identihed by eye 
and measured by contours down to a limiting sur- 
face brightness about lOcr above the background. 
We single out 30 star-forming complexes in this 
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Fig. 8. — i? 435 — Vgoe vs. Veoe— Fsi 4 colors of all SCCs in HCG 59, plotted along with i Zq (solid red line, 
bottom) and Zq (dashed red) model tracks by Marigo et al. (2008), spanning an age range of 6 Myr to 
13 Gyr (indicated by filled stars). The green lines (solid and dashed for | Zq and Zq) show a Starburst99 
track that accounts for nebular emission, as suited to the very youngest clusters that have not yet expelled 
their natal gas. Data-points that lie to the left of the dotted green line are considered ‘nebular’, i. e. younger 
than ^ 10 Myr. The high-confidence sample of SGGs (those with My < —9) is denoted by solid dots, 
while open triangles show fainter candidates comprising the extended sample (all cluster candidates to the 
detection limit). Much of the spread in the extended sample colors is due to stochasticity in populating the 
stellar initial mass function within clusters, as discussed in the text. We also indicate an extinction vector 
of length Aqqq = 1 mag. We will present the distribution of each galaxy separately in the next figure. 
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Fig. 9 . — S435— Veoe Vgoe— -^8i4 colors of SCCs, divided by galaxy, following Figure 8 . Typical error- 
bars are indicated as crosses in the lower left part of each plot. The irregular galaxy D is the only one 
that displays pronounced current star formation activity, as deduced by the number of nebular sources - 
although more nebular sources could be hidden behind dust and gas in the highly inclined spiral galaxy C. 
Galaxies C and D are consistent with continuous star formation histories over the course of their histories. 
The absence of GCs in D might indicate the galaxy is young, although it could be due to the fading of its 
low-mass cluster population below our detection limit (which is a function of age). Galaxies A and B only 
appear to host globular clusters and a handful of intermediate-age and young clusters, indicating quiet star 
formation histories over the past Gyr or so. 
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Fig. 10. — Veo6 vs. Veoe— Fsi 4 color-magnitude diagrams, following the datapoint plotting conventions of 
Figure 8. We have provided Marigo et al. (2008) model tracks for four different cluster masses, times two 
metallicities, | Zq and Zq (top and bottom) (solid lines), with the ‘nebular’ tracks plotted as dotted lines. 
As in Figure 9, the age increases from top left to bottom right. In both samples, we tentatively deduce typical 
cluster masses - with the caveat that mass and extinction are degenerate - of < 10^ Mq, with an additional 
focus of possible 10® Mq clusters. This distribution is consistent with local cluster populations, given the 
imposed cut-off for the high-confidence sample at My < —9 (and therefore the lower limit of ^ 10^ Mq). 
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Fig. 11. — Veoe 'vs. V 606 ~^ 8 i 4 color-magnitude diagrams for each galaxy individually, following the plotting 
conventions of Figure 10. Galaxy B is notable for the large population of apparently old, massive (~ 10® Mq) 
clusters, while galaxy D has shows a sizable population of 10^ Mq clusters consistent with a predominantly 
young and intermediate-aged system. 
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compact group, three in galaxy C, and 27 in D. 
We find no complexes in the E/SO galaxy B or A, 
the evolved (Sa) spiral. The sizes of these com- 
plexes are comparable to those in local systems 
(Figure 12; cf. Elmegreen et al. 1994, 1996) and 
consistent with HCG 7, a relatively inactive com- 
pact group (Konstantopoulos et al. 2010). Their 
colors are indicative of star formation and nebu- 
lar emission, as expected of star-forming regions 
(Figure 12, right). In all, the star-forming com- 
plexes in HCG 59 are comparable to their counter- 
parts in local star-forming galaxies. They follow 
very closely a relation between brightness and 
size, much more so than we found in our previ- 
ous study of HCG 7. We show this correlation in 
Figure 12 (right) where the line is a simple linear 
fit to the data. The linear fit implies a similar 
surface brightness for complexes across the group, 
unlike in HCG 7, where we found significantly 
more scatter. 

3.3. The ancient globular cluster systems 
of HCG 59 

Globular clusters represent the earliest eras of 
star formation in a given galaxy. Their color dis- 
tributions carry the imprints of interactions and 
mergers and thus may probe the history of their 
hosts over a very long timescale. 

In this section we provide a full analysis of the 
number and color distributions of the GC popula- 
tions in HCG 59 and complement the star forma- 
tion and interaction histories we began to explore 
in Section 3.1. This analysis is built around the 
plots of Figure 13. The top row has a color-color 
plot (left) and histogram of the {B — I)q distribu- 
tions (right) in each galaxy. (The color-magnitude 
diagram is shown in the left panel Figure 3.) Ow- 
ing to their projected proximity, the GC systems 
of galaxies A and D overlap, and so we cannot pro- 
vide separate analyses. However, given the rela- 
tive masses of the two galaxies and their evolution- 
ary stages, it is safe to assume that the population 
is dominated by A. 

We find a fairly large population of globular- 
clusters across this compact group, the vast major- 
ity of which are found in and around B. Recall that 
given the fading of clusters with age, GCs need 
extremely high masses, M > 10® Mq, to be de- 
tected to such large distances. We also find some 
clusters at large radial distances from the center of 


this galaxy, including some found along the stellar 
stream that appears to connect galaxies A and B 
(see Section 4.3) and its projection on the far side 
of galaxy B. Galaxies A and C host small pop- 
ulations. To quantify these populations, we de- 
rive the specific frequency, 5'at, a measurement of 
the number of clusters per unit galaxy luminos- 
ity, for each galaxy. First, we correct the observed 
number of GC candidates by the background cor- 
rection noted above, and then calculate the total 
number of GCs expected around each galaxy by 
first adopting a photometric completeness fraction 
of / = 0.9 ± 0.1 for objects with Veoe < 26, and 
correcting for the expected fraction of GCs that lie 
below this magnitude limit. The Sm values are 0.3, 
7.7 and 0.1 for galaxies A, B and C respectively, 
assuming that all bright GCCs consistent with the 
halo of A are actually bound to A. The measured 
and derived numbers are collected in Table 4. 

The size of the population in galaxy C is consis- 
tent with its Sc morphological type, while B has a 
tremendously rich system, about twice the number 
of GCCs expected. Conversely, the GC population 
of galaxy A is much poorer than expected, com- 
pared to the values of Sn 1 typically seen in 
Sa spirals (Chandar et al. 2004). Many of the de- 
tected GCCs appear to lie in a ring just outside of 
the bulge of the galaxy. Comparing to other GC 
populations in HCGs, our specific frequency for 
HCG 59B of S' AT = 7.7±3.0 is larger than (yet still 
consistent within the uncertainties) those observed 
so far in other large elliptical galaxies in HCGs: 
global values of Sn = 3.6 ±1.8 for HGG 22A from 
Da Rocha et al. (2002),^ and Sn = 4.4 ± 1.3 for 
HGG 90G from Barkhouse et al. (2001). 

Inspection of the region marked ‘outer B’ in 
Figure 3 shows a possible excess of GC candi- 
dates over the background level described above 
in Section 2.4, perhaps tracing an intra-group stel- 
lar population. This is not unusual, as compact 
groups by their very nature are likely to promote 
interactions. (White et al. 2003) detected diffuse 
intragroup emission in HCG 90 accounting for up 
to half the light of the group. Two further stud- 
ies, Da Rocha & Mendes de Oliveira (2005) and 
Da Rocha et al. (2008), found intra- group light in 


'^The lower Sjv ~ 1-7 ± 1.0 from Barkhouse et al. (2001) for 
HGC 22 A is due to their larger My for the galaxy - both 
studies find a very similar total number of globular clusters. 
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Fig. 12. — Color-color and size-luminosity {D — Mb) diagrams for star cluster complexes in HCG 59. The 
plotting of the evolutionary tracks follows the conventions of Figure 8. Complex symbols are coded to 
indicate where they are found: galaxy C (filled, yellow triangles) or D (filled, orange circles) The colors 
of the complexes are indicative of star formation, as expected. We find the complexes to follow a linear 
luminosity-size relation (see the upper left corner of the right-hand panel) very well. This implies a similar 
surface brightness for complexes across the group, unlike in HCG 7, where this relation is not obeyed. 



Fig. 13. — Left: Color-color diagram of all detected GCCs. The dotted parallelogram on the left shows the 
selection box, based on the colors of GCs in the Milky Way (Harris 1996). Center: {B — I)q distribution 
for all detected GCCs. The shaded area in the plot of the A/D distribution omits GCCs superposed over 
the body of D. The ‘outer B’ distribution covers all GCCs in pointing 2 that are at least 15 kpc from the 
centre of galaxy B. This region might suffer from some contamination, although it is consistent with observed 
distributions of GCs, which show a tendency for red clusters to be at large galactocentric distance as would 
be expected for a halo population. On the right, we relate the color distribution to that of metallicity and 
compare galaxies A and B directly. The dotted line in the GC metallicity distribution of galaxy B shows a 
best fit Gaussian, which we plot since a more complex distribution is not statistically justified. 
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another six groups: HCGs 15, 35, 51, 79, 88, and 
95. In HCG 59, a total of iV = 15 ± 4 GGGs lie 
within the 2.7 arcmin^ region, where we expect a 
background contribution of 6 ± 2 objects. Assum- 
ing only Poisson noise, this suggests a ~ 2 (t excess 
of objects in this part of HGG 59. This possible 
excess could be due to variations in the stellar den- 
sity in the Galactic halo or the Sgr Stream. How- 
ever, the luminosity function of this small number 
of sources in the ‘outer B’ region (shown as the 
solid line in the right panel of Fig. 3) is weighted 
towards the faint end, suggesting that at least 
some of these objects are indeed true GGGs lo- 
cated far 25 — 50 kpc) from galaxy B. The un- 
usually large population of GGGs associated with 
B, the anomalously small population in A, and the 
possible population of GGs in the IGM (presum- 
ably stripped from a member galaxy) , all indicate 
that this compact group environment may have re- 
distributed GGs between member galaxies and/or 
to the IGM. 

In particular, possible interactions in the re- 
cent history of galaxy B might have redistributed 
its GG system. With this in mind, we compare 
the azimuthally averaged radial profile of the spa- 
tial distribution of the GG system with the surface 
brightness profile of the galaxy, as derived through 
GALFIT (Peng et al. 2010b). The galaxy is best 
fit by a single-component Sersic profile of shape 
n = 3.1. Interestingly, we find the Sersic profile of 
the galaxy to provide a better description of the 
GG system than the best fitting power-law profile 
of index —1.3. This is contrary to the finding of 
power-law GG system distributions around loose 
group member NGG 6868 and HGG 22 A. If that 
is to be considered the norm for compact groups, 
then perhaps the recent interaction activity about 
galaxy B has changed the shape of its GG halo. 

The bottom panel of Figure 13 shows the {B — 
I)o color distribution (converted to the John- 
son photometric system for direct comparison to 
Galactic globulars) of the clusters in each galaxy. 
The shaded area in the plot of A/D shows clusters 
that are clearly part of A, i. e. omits the ones that 
are projected upon the body of D. This does not 
alter the distribution, strengthening our assump- 
tion of a small population in D. Galaxy G also 
hosts a very small population, as expected due to 
its low mass. 

The color distribution can act as a proxy of 


metallicity for GGs and we take advantage of 
that to compare the two populous distributions 
of galaxies A and B. The color distribution of 
galaxy A seems fairly flat and the low numbers do 
not allow for a statistical treatment. Galaxy B, 
however, provides a large enough population to 
perform a test for bimodality, using the KMM al- 
gorithm of Ashman et al. (1994). This returns no 
evidence for a composite distribution in the metal- 
licity distribution. 

4. Discussion 

4.1. Extending the membership of HCG 59 

In Section 2.5 we introduced a search for dwarf 
galaxies in HGG 59, which we continue here. Since 
all objects we are considering here are covered 
by SDSS, we will not provide images and spec- 
tra here. More information can be obtained from 
the SDSS database using the plate IDs, Modified 
Julian Dates (MJDs), and fiber IDs given in Ta- 
ble 3. 

Regarding the morphologies of the new mem- 
bers, I, the candidate covered by our HST imag- 
ing, seems irregular, with a peaked light profile. 
This agrees with its spectrum, which shows clear 
emission lines and a continuum shape typical of a 
spiral. We classify it as dim. The rest of the galax- 
ies are not covered with high-resolution imaging, 
so we are more conservative in classifying them. 
F shows quite clear spiral structure in the SDSS 
images, and is therefore given an Sd type. We 
note that its star formation seems to be declin- 
ing, given that Ha is the only detectable emission 
line. G and H appear quite irregular and elon- 
gated, with spectra exhibiting blended emission 
and absorption. We assign them dir types. Fi- 
nally, J shows a morphology closer to spherical 
and weak emission lines (although the S/N does 
not allow for a confident determination); we as- 
sign it a dE type. 

Regarding their roles as group members, four 
of the five galaxies do not appear to be interact- 
ing with any other members, as might be expected 
from their locations far from the group core (Fig- 
ure 14, left). It is only I that shows some evidence 
of an interaction with galaxy B, in the form of the 
‘B-I arc’ described in the Section 2.5. This stellar 
stream includes the large star-forming region we 
find to the west of galaxy B, which is part of the 
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Table 4: Observed and derived properties of the globular cluster systems in HCG 59 


Galaxy 

Nggc^ 

A^back'^ 

Nv<26^ 

A^total'" 

Sn ^ 

A 

18±4 

6±1.0 

13±6 

33±21 

0.3T0.2 

B 

162±13 

6±1.0 

180±34 

462A182 

7.7T3.0 

G 

7±3 

6±1.0 

1±3 

3±8 

O.liO.3 

D 

11±3 

1±0.2 

11±5 

28±18 

1.6T1.0 


“Number of detected GCCs 
^'Background correction 

“Number of clusters brighter than V = 26 mag 
“^Estimated total GC population 
“Specific frequency 


analysis of Section 3.1. 

In order to assess whether these galaxies belong 
to the group, we perform a phase-space analysis, 
following the statistical studies of Mulchaey & Zabludoff 
(1998) and Zabludoff & Mulchaey (2000). This 
is shown in Figure 14, right panel, where each 
datapoint represents a dwarf galaxy. The x-axis 
measures the distance of a galaxy from the group 
centroid, calculated as the mass-weighted average 
position of the four main members (A through 
D). The y-axis shows the offset of a galaxy’s ra- 
dial velocity from the group mean, normalized to 
the core group (A through D) velocity dispersion 
of 314 kms“^. 

We find all galaxies to satisfy our member- 
ship constraints. Four of five spectroscopically 
detected galaxies lie within the boundary set by 
the four main members; galaxy G, which has the 
largest offsets in physical and velocity space, is 
moving with a radial velocity offset less than three 
times the group dispersion. We therefore con- 
sider all five galaxies under consideration here to 
be members of HGG 59, based on a strict 3cr- 
clip, as demonstrated in Figure 14. The inclu- 
sion of new members updates the velocity disper- 
sion of the group to 335 kms“^. Based on this 
value, we derive a dynamical mass for HGG 59 
of Mdyn = 2.8 X 10^^ Mq, a ~ 10-fold increase 
with respect to the main members alone. This 
in fact changes the J07 evolutionary stage of the 
group from Type II (intermediate) to Type III 
(gas-poor), as it yields a ratio of H i-to-dynamical 
mass of 0.71, with the caveat that the measured H i 
mass is likely underestimated because the dwarfs 
at large group radii are not included. 


It is interesting to find dwarf galaxies at large 
distances from the center of the group. This lack 
of barycentric clustering is also observed in the 
Local Group, where it is seen as a morphology- 
density relation: dwarf irregulars (dir) are found 
at larger distances from the group center (Grebel 
1999) than the quiescent dwarf spheroidals (dSph) 
and dwarf ellipticals (dE). This may indicate that 
some dirs are galaxies experiencing their first in- 
falls to the group center. Such a situation could 
explain the relatively large velocity offset and the 
star- forming nature of galaxy G. 

The membership of galaxy G seems the most 
uncertain of the hve galaxies discussed above, 
given the marginal agreement with the 3 a velocity 
cut, and the large projected barycentric distance. 

This is important, as its inclusion does affect the 
updated dynamical properties, due to the large 
change in group radius. We quantify this in Ta- 
ble 5, where we summarize the dynamical proper- 
ties of HGG 59. Those numbers show a change of 
mass by a factor of ~ 4 or 10 with and without 
galaxy G, while the J07 type changes from II to HI 
regardless of the inclusion of G. The velocity dis- 
persion is most affected: including G increases the 
value to 336kms“^ from the original 314kms“^, 
while excluding G significantly reduces the disper- 
sion to 208kms“^, which is more consistent with 
the evolved state suggested by our analysis. 

4.2. The current state of star formation 

We have presented several diagnostics of star 
formation activity across HGG 59. The Tzanavaris et al. 
(2010) SFRs of galaxies A through D are ~ 
[4.99,0.02,0.16,0.48] Mq yr“^. These are de- 
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Fig. 14. — Diagnostic diagrams of dwarf galaxy membership. The left panel shows the positions of the dwarfs 
and giant galaxies as blue dots and red stars respectively. We have labelled the dots with the appropriate 
lettering for HCG members. On the right we present a phase-space diagram that plots the offset of a 
galaxy’s radial velocity from the established group mean (cTcore, the line of sight velocity dispersion of 
the giant galaxies) against its distance from the group barycenter (mass- weighted centroid). Four of five 
spectroscopically detected dwarfs lie within the boundaries set by the main galaxies and as such are confirmed 
as members. The fifth is moving at a relative radial velocity less than three times the group’s velocity 
dispersion and is therefore also included as a member. The lower boundary of this 3 cr-clipping criterion is 
shown on the plot as the shaded area. 


Table 5: Updated dynamical properties of HCG 59 


Reference 

V 

( kms 

(7y 

^^dyn 

(Mo) 

J07 type 

Hickson (1982) 

4036 

314 

2.9 X IQi^ 

II 

This work 

3924 

336 

2.8 X 1Q13 

HI 

This work, excluding G 

4015 

208 

1.1 X 10^3 

HI 
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termined from the combination of UV and IR 
light assuming that all of the light emanates from 
star-forming regions. The presence of young star 
clusters in galaxies C and D shows that stars are 
forming at a fair pace, in accord with the star 
formation rates quoted in Table 1. Overall, the 
SFRs are consistent with the infrared SEDs of 
these galaxies with the notable exception of A 
which is almost certainly strongly contaminated 
by AGN emission and shows no evidence of on- 
going star formation from the other evidence on 
hand. There are several soft X-ray point sources 
throughout the group, which are likely to probe 
compact stellar remnants local to HCG 59 galax- 
ies. There is also soft, diffuse X-ray emission con- 
fined to the galaxies, which probes ~ 10® K gas 
heated by star formation (stellar winds and SNe). 
The IR images do not show much that is surpris- 
ing: emission along the spiral arms of C and in 
the star-forming clumps embedded in D. Four of 
the five dwarf galaxies show star formation activ- 
ity, J being the exception. They are found to be 
star-forming based on either their bright emission 
lines or blue continua. Galaxies A and B, on the 
other hand, exhibit quiescent or even extinguished 
star formation. In the case of A, this is inferred 
by the absence of young star clusters. 

In all, the group does not appear to be un- 
dergoing a burst or any other event notable in 
terms of current star formation. With the ex- 
ception of the two most massive galaxies, the 
group is forming stars at a regular pace. This 
is also exemplified by the behavior of star clus- 
ter complexes across the group, which follow very 
closely after a brightness-radius relationship con- 
sistent with typical nearby galaxies, and in con- 
trast to HCGs 7 and 31 (Konstantopoulos et al. 
2010; Gallagher et al. 2010). 

4.3. Signatures of interactions in the intra- 
group medium 

Major interaction or merger events are very 
often accompanied by bursts of star and cluster 
formation. The examination of star clusters in 
HGG 59 presented in Section 3.1 did not show any 
evidence of such events in the last few Gyr. Given 
the high mass-detection limit for star clusters at 
this distance, they cannot be used to trace minor 
dynamical events. We therefore search for such 
evidence in the lowest surface brightness features 


detectable in our LCO images. In Section 2.5, we 
reported the detection of a low surface brightness 
stream of material in the projected area between 
galaxies A and B. This is visible at low surface 
brightness in our HST imaging and quite pro- 
nounced in the wide-field images from Las Cam- 
panas, at the > 3 ct level. We also detected an arc 
of luminous material to the west of galaxy B, per- 
haps connecting it to compact galaxy I, which we 
discuss in Section 4.1. 

The low surface brightness and limited extent 
of the ‘B-I arc’ preclude precision photometry. We 
can therefore only pursue an in-depth analysis of 
the bridge between A and B. In order to derive 
the photometric properties of this feature, we first 
drew a color-map to look for an evident B — R 
gradient. Unfortunately, the stream is not bright 
enough to clearly dominate the image background. 
We thus conducted photometry of the area and 
the outskirts of the two galaxy that bracket it. 
We used large apertures of radius 75 pc (8 px) 
to reduce the background noise. The results are 
plotted in Figure 15; on the left panel we compare 
the measured photometry with the evolution of 
the B — R color, according to the Marigo et al. 
(2008) model tracks of three metallicities. We find 
an intermediate value between the colors of the 
outskirts of the two galaxies that define this region 
(plotted as yellow dashed lines), cautioning that 
the emission in this region might be affected by 
the two galaxy light-envelopes to some extent. 

The origin of the bridge is not clear and its 
faintness makes it difficult to ascertain the domi- 
nant source of emission. It could consist either of 
stripped stars, or stars that formed in situ from 
stripped gas. If this is mixed stellar material from 
the two galaxies, we cannot study it in any more 
detail. We can, however, develop the in situ for- 
mation scenario further, by treating the bridge 
as a simple stellar population. In this case, the 
color-magnitude diagram of Figure 15 (right), pro- 
vides an age estimate of about 1 Gyr, depending 
on metallicity. In addition, the GMD plotted in 
this figure provides an estimate of the stellar mass 
contained: with 10^ Mq of ~ 1 Gyr old stars in 
each aperture, we extrapolate a mass in the order 
of ~ 10® Mq, i. e. a density of 100 Mq pc^. 
If the stars here are stripped from a galaxy, the 
overall mass will be higher, as the M/L of simple 
stellar populations increases with time. 
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Fig. 15. — Diagnostic plots of the tidal bridge between galaxies A and B. The age-color plot on the left 
shows the evolution of the B — R index over time (logr), for Marigo et al. (2008) models of three different 
metallicities: | and solar, shown as dashed green, dashed-dotted blue and solid red lines respectively. The 

range of colors in the bridge, as measured in five successive apertures and shown as a shaded region, lines 
up with the tracks at ages around 1 Gyr (depending on metallicity). The photometry of this faint feature 
might be affected by diffuse light from the galaxies, which brackets that of the bridge (denoted by dashed 
yellow lines, galaxy A is the bottom line). A secondary diagnostic is presented on the right, in the form of 
a color-magnitude diagram. The numbered dots denote log age in years, while the stars show the color and 
magnitude of each bridge aperture. If this is a simple stellar population (see text), then we can extrapolate 
a mass of ~ 10® Mq over the 1 kpc^ area of the bridge. This measurement can provide constraints for 
models of the evolution of tidally induced stellar structures. 


28 


4.4. HCG 59B as a merger remnant 

HCG 59B, the E/SO galaxy on the west side 
of the group, seems quite regular at first glance, 
however, a close inspection of the low-level light re- 
veals some interesting features. While ELLIPSE 
htting shows an overall smooth isophotal struc- 
ture in B435 light, there is severe isophotal twist- 
ing in the central regions. This is spatially co- 
incident with several patches of extinction we de- 
tect in the HST images; they are most pronounced 
in i?435, observable in Veoe and hardly detectable 
in /si4, implying a thin column of dust. In the 
Spitzer bands, the fits indicate very symmetric 
structure in the 3.6 and 4.5 fim bands, but the 5.8 
and 8.0 /rm fits show a faint cusp some ^ 7 px, or 
^ 2.5 kpc from the surface brightness profile 
peak. 

The color-composite IRAC image (Figure 5) 
shows hints of structure in galaxy B; however, 
an evolved elliptical/lenticular with near-zero star 
formation should present a smooth isophotal pro- 
hle in all bands. Nonetheless, its SED (Figure 6) 
shows a gradual decline, indicative of emission 
from stellar photospheres, rather than the heated 
dust associated with star formation. 

To investigate these irregularities further, we 
take advantage of the high resolution of the HST 
images. We construct pixel-by-pixel color maps 
of B, with the aim of tracing the exact location 
of the extinction patches. If the underlying stel- 
lar population is evolved to the same stage {i. e. 
an evolved SSP), then extinction will be the only 
source of discrepancies in color. There are, in fact, 
three possible sources of color variations in this fil- 
ter combination: (i) extinction across a similarly 
colored stellar population (mixed or coeval); (ii) 
spatially separated stellar populations of various 
ages and/or metallicities; or (iii) the presence of 
gas and star formation - i. e. H ii region emission 
lines. 

The maps, shown in Figure 16, cover the three 
possible filter combinations; for reference, the Veoe 
image is also shown. The B — I map is the 
one most sensitive to extinction. To quantify, 
the Cardelli et al. (1989) extinction law assigns al- 
most twice the extinction in the B435 as it does 
to the /si4, A435/A814 = 1.85 (A43s/Ay = 1.13, 
cf. Asia! A y = 0.61). This is therefore the map 
we use to detect patches of low extinction, of 


order 0.3 mag, in three hngers extending ap- 
proximately eastwards from the north-south line 
through the nucleus. This faint structure is seen 
in all three colormaps, but not in the Veoe image. 
The F606W filter covers various emission lines, in- 
cluding Ha, [Nil], [Sil], H/3, and [Oiiij. Interest- 
ingly, Martinez et al. (2010) found Ha, [Nil] and 
[SH] emission in the spectrum of B, at relative in- 
tensities consistent with a composite H ii region 
plus AGN emission. From the concentrated, blue 
core of the V — I image, this line emission appears 
to be spatially coincident with the nucleus. 

Gircumstantial evidence for a close encounter 
in the recent history of galaxy B is provided by 
the uneven distribution of its GG system. GGs 
normally form spherical haloes, however, here we 
hnd GGs at large radii, many concentrated along 
the stellar stream that seemingly connects galax- 
ies A and B, and its extension across the far side 
of galaxy B. Furthermore, the overabundance of 
GGs in galaxy B is matched by a severe dearth 
of clusters in A. Given the possibility that the 
two galaxies interacted ~ 1 Gyr ago, a scenario 
whereby GGs are transferred between the two sys- 
tems is not out of the question. It is unclear from 
a dynamical perspective why in the process of an 
interaction the GGs would flow from A, the more 
massive entity, to B. In a simple thought experi- 
ment, we move as many clusters from B to A as are 
required to level the Sn of A to the nominal value 
for an Sa. This still leaves an excess of GGGs in 
B relative to normal. However, the factor-of-two 
uncertainties involved in the determination of Sn 
do not rule out this scenario of GG ‘swapping’. 

The X-ray map of this galaxy, as described in 
Section 2.7, reveals two distinct X-ray sources in 
the nuclear region. Unfortunately, with full-band 
luminosities of Lx = (1.4, 1.7) x 10^® erg s“^, nei- 
ther is sufficiently luminous be identified as an 
unambiguous AGN - a possibility that the opti- 
cal spectroscopy of Martinez et al. has indicated. 
Given the lack of ongoing star formation in the 
region, the sources are unlikely to be high mass 
X-ray binaries, though luminous low mass X-ray 
binaries (associated with older stellar populations) 
or groups of them unresolved at the distance of 
HGG 59 are plausible. In addition, due to the 
uncertainty in matching X-ray sources to optical 
imaging, we cannot confidently derive a one-to-one 
correlation between the optical clumps and these 
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Fig. 16 . — HST color maps of the central region of HCG 59 B. The maps are centered on the nucleus, defined 
as the peak of the /si4 surface brightness profile, and aligned to the world coordinate system, with north 
to the top and east to the left. The bottom right map shows the Veoe brightness. Contours show .B435, 
^606; -^814; ^nd /gi4 flux in the four panels (upper left, upper right, lower left, lower right). Each successive 
contour maps a difference in brightness of 0.5 mag. The B — I map on the bottom left is most sensitive to 
extinction and does indeed show patches of moderate reddening. The top row maps also show the extinction 
patches evident in the bottom row. In both the V — I and B — I color maps, the nucleus appears bluer than 
the surroundings; this is consistent with additional flux from line emission in the .B435 (H/ 3 ) and Vgoe (Ha) 
filters. The meaning of the observed colors is discussed in more detail in Section 4 . 4 . 
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sources, although the correlation is confirmed to 
within the X-ray positional uncertainties. We note 
however that, as discussed in Section 2.7, the nu- 
clear X-ray sources are spatially close to the opti- 
cal center of the galaxy. 

The excess of globular clusters and uneven dust 
distribution in the nuclear regions of B hint at 
some interaction in the more distant (> Gyr) 
past, but the lack of additional evidence for struc- 
tural disturbances limits our ability to infer more. 
We do detect a few young clusters in galaxy B 
and a non-zero (though low) SFR, and therefore 
some reservoir of cold gas is present. Accretion 
of a satellite galaxy is therefore a possibility. 
Furthermore, the unimodal GC color distribu- 
tion does not favor a gas-rich, major merger in 
the past. This conclusion follows the paradigm 
of Muratov & Gnedin (2010), who attribute the 
known bimodality of GC colors to late-epoch 
mergers. It is also a reasonable assumption that 
the implied interaction did not feature a major 
merger with a gas-rich system, as that would have 
enhanced the young and intermediate-age cluster 
populations. 

4.5. Nuclear activity in HCG 59A 

In Section 2.6 we reported that the IR emission 
in galaxy A is consistent with being dominated by 
an AGN, rather than star formation. This is based 
on the disparity between the galaxy’s morphologi- 
cal type of Sa and the high UV-I-IR SFR suggested 
by interpreting the emission as related to SF. Fur- 
thermore, the lack of young massive clusters is in- 
consistent with a SFR of 5 M0yr“^. 

This hypothesis is supported further by the 
hnding of a hard X-ray source in the nuclear re- 
gion of galaxy A, with Lx = 1.1 x 10^° erg s“^ as 
reported in Section 2.7. The spectroscopic AGN 
survey of HCGs (Martinez et al. 2010) places the 
galaxy at the interface of the H ii and AGN zones 
in the ‘BPT’ diagram (Baldwin et al. 1981), based 
on optical emission-line ratios. The 2"-wide slit 
they used encompasses a wide region (effective 
aperture of 0.58 kpc), therefore the signal is most 
likely diluted by circumnuclear and disk light. 
Visual inspection of the central region reveals 
asymmetric structure, resembling a second nuclear 
source of comparable I8i4 luminosity to the nu- 
cleus. As in the previous section, we employ HST 
color maps to take advantage of the spatial reso- 


lution of ~ 12 pc per pixel. 

The maps, shown in Figure 17, reveal a cone 
of blue light at {x,y) ~ (—0.1,0.!), with colors 
of i?435 — Vqoq ~ 0.9, Veoe ~ .^8i4 — 0.6, and 
B 435 — Isi 4 — 1-6. The complexity of the central 
region inhibits easy interpretation, but one clue 
are the emission lines covered in the three bands. 
Veoe covers H/3, [O iii], Ha and [N ii], which can be 
associated with star formation and/or AGN activ- 
ity. In this scenario, the blue cone could stand out 
as a result of geometry, perhaps being located in 
a break in the dust distribution. A simple expla- 
nation could relate this feature to an unreddened 
line of sight through the inner spiral structure of 
A. In that case, however, we would expect to see 
bright, young star clusters, as they are known to 
shine through thick columns of dust, let alone rel- 
atively dust-free regions (e. g. Region B in M82; 
Smith et al. 2007; Konstantopoulos et al. 2008). 

A different interpretation can relate this struc- 
ture to an AGN. The lack of symmetry could 
suggest a small narrow-line region, photoion- 
ized by the AGN continuum, with projection ef- 
fects and obscuration hiding the cone on the far 
side. This would produce strong [O ill] and [N ii] 
emission, the presence of which was reported by 
Martinez et al. (2010). This geometry is consis- 
tent with the inclination of galaxy A of no more 
than 30° (assuming the AGN and galaxy share 
the same inclination angle). 

Combining the pieces of evidence collected from 
the X-ray, optical and MIR emission, we propose 
that the nuclear emission in HCG 59A is domi- 
nated by a low-luminosity AGN with a photoion- 
ized narrow-line region. The onset of activity may 
be related to a possible encounter with galaxy B 
about 1 Gyr ago, as tentatively dated from the 
colors of the bridge connecting the two galaxies 
(Section 4.3). 

5. Summary 

We have presented an analysis of HCG 59, 
a compact group comprising four main galax- 
ies and at least five newly discovered dwarfs at 
the Mr < —15.0 mag level. Our results are 
based on multi-wavelength observations and con- 
tinue a series of papers that have followed two 
different approaches: on the one hand we have 
treated the overall properties of Hickson compact 
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Fig. 17. — HST color maps of the nuclear region of HCG 59A, following the plotting technique presented in 
Figure 16. North is to the top and east is to the left. The central peak appears elongated in the Veoe 
map (lower right) and shows a composite structure with two additional peaks nearby and a long tail. The 
tail could be an artifact introduced by the complex and heavy extinction in this region, perhaps caused by 
the inner spiral structure of HCG 59A. It could alternatively be associated with the nuclear activity we note 
in the X-ray image of this galaxy (Section 2.7), as well as the significant excess in 8 and 24 ^m emission, 
discussed in Section 2.6. A blue cone-shaped structure is most evident in the V 606 ~-^ 8 i 4 color map (top right). 
The nature of this blue cone is discussed in Section 4.5. 
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groups (Johnson et al. 2007; Gallagher et al. 2008; 
Tzanavaris et al. 2010; Walker et al. 2010); on the 
other hand, we have surrounded our HST obser- 
vations with a multi-wavelength dataset to pursue 
in-depth investigations of individual CGs, one at 
a time (Palma et al. 2002; Gallagher et al. 2010; 
Konstantopoulos et al. 2010). 

Gompared to HCGs 7 and 31, two compact 
groups previously studied in this series, HGG 59 
presents something of an intermediate step: where 
HGG 7 was found to be interacting solely in the 
dynamical sense {i. e. currently in the absence of 
direct hit encounters), HGG 59 shows evidence for 
stronger interactions in the recent past. There 
is evidence for star formation in the intragroup 
medium in the H ii regions to the Northwest of 
galaxy B, in sharp contrast to HGG 7 where the 
star-formation associated with each galaxy was 
self-contained. In the context of the evolutionary 
sequence we proposed in Konstantopoulos et al. 
(2010), it occupies a stage further along than 
HGG 7. It has begun building an IGM, as tes- 
tified by some amount of intra-group light, but 
has yet to build up a large elliptical fraction. Its 
classification as a relatively unevolved group (J07) 
is in accord with standard diagnostics such as its 
velocity dispersion of ~ 335 kms“^ and lack of 
diffuse, extended X-ray emission. 

Through the use of SDSS data, we associated 
five dwarf galaxies to HGG 59 for the first time. 
Their inclusion updated the velocity dispersion 
and dynamical mass of the group and changed 
its JO 7 evolutionary type to an evolved group of 
Type HI (originally Type II, or intermediate). The 
star-forming nature of these dwarfs, their radial 
velocities, and distances from the group core seem 
to suggest that some may be infalling for the first 
time. The noted relation between morphological 
type and barycentric distance follows the one ob- 
served in the Local Group and may be considered 
- conversely - as an indication of the ‘young’ dy- 
namical state of HGG 59. The above information 
highlights the importance of studying the dwarf 
galaxy contingent of compact groups. 

Star formation is proceeding at a regular pace 
in this GG, certainly at a rate consistent with the 
morphologies of the member galaxies. The reg- 
ularity of star cluster complexes agrees with this 
image. The star cluster population does not show 
evidence of major, gas-rich, interactions in the 


past Gyr and the IR SEDs are generally as ex- 
pected. One exception to this rule is galaxy D. 
Given its large size, the morphological regularity 
of its neighbors, and the lack of evidence of recent 
interactions, it is not clear why it has such an ir- 
regular structure and why it is forming stars at 
the rate that it is. 

Where the information is unclear for galaxy D, 
B is evidently in the midst of at least one dy- 
namical process. This probably started more than 
a Gyr ago given the lack of tidal features such 
as shells and tails commonly observed in such 
events in the optical (Schweizer & Seitzer 1998; 
Hibbard et al. 2001; Mullan et al. 2011). It is 
likely physically associated with an arc of star 
formation that seemingly connects it to dwarf 
galaxy I (the arc that hosts the discovered ex- 
tragalactic Hii regions), and 59B lies at one end 
of a low surface brightness bridge of stellar ma- 
terial, which might physically connect it with A. 
We found this feature to be no older than ~ 1 Gyr 
by age-dating its stellar population. Perhaps most 
intriguingly, the globular cluster population of B 
is anomalously large, with significantly more glob- 
ular cluster candidates than A, despite its lower 
stellar mass. The origin of this discrepancy is un- 
clear, but hints at an additional event in the more 
distant past. 

In one sense, this group is typical of early or 
intermediate-stage HGGs. There are plenty of dy- 
namical processes at play, however they are all 
proceeding at a low level and centered around one 
or two objects. Galaxy B is at the focus of all 
such processes that our diagnostics can reach, like 
galaxy B in HGG 7 (also an early-type galaxy). 
In addition, both these groups may feature an ac- 
tive nucleus host galaxy (the dominant galaxy A 
in both groups, albeit the detection is not certain 
in HGG 7). 

We find the presence of a low-luminosity AGN 
in A to be likely given the spatial coincidence of 
the galaxy centroid with a ~ 10^° erg s“^ X-ray 
point source. The inferred SFR from the IR-I-UV 
luminosity is not supported given the lack of young 
star clusters, and therefore likely is overestimated 
because of significant AGN contamination. If the 
A-B bridge constitutes a physical connection be- 
tween the two galaxies, a causal connection be- 
tween the interaction and the AGN is possible. In 
this case, we can constrain the timescale since the 
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onset to no more than one Gyr. 

In the introduction, we discussed HCGs as a 
potentially special environment in terms of galaxy 
evolution. The apparent duality of ‘modes’ in 
which HCG galaxies are found - either star- 
forming or quiescent - and the evident lack of an 
intermediate stage population are in accord with 
the mid-IR color-space ‘gap’ discussed in our pre- 
vious work (Johnson et al. 2007; Tzanavaris et al. 
2010; Walker et al. 2010). The impact of the 
group environment will be the topic of the 
next paper in this series. There we will treat 
HGGs 16, 22 and 42 with a goal of understand- 
ing the evolutionary processes at play in compact 
groups and relating their galaxies to those found 
at other levels of galaxy clustering. 
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